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The design of a flight vehicle is a very complex problem. Technical, 
economic, marketing and other challenges must be overcome to achieve success. 

It is not surprising that "design 1 ' means different things to different people. 

To the aerodynamicist it usually means "given a set of criteria, find the best 
configuration to achieve them." Even if we limit our consideration to "aero- 
dynamic design," we are still faced with a complex problem. A team of people 
must use a combination of intuition, experience, testing, numerical and 
theoretical analysis to develop a successful design. The numerical part of 
this process has improved substantially in the past decade due to the develop- 
ment of computational methods which can treat realistic geometries and 
governing equations. A perspective of this development is given by Graves . 

This paper discusses one computational approach to aerodynamic design, namely 
the procedure of using numerical optimization methods coupled with computational 
fluid dynamic (CFD) codes. 

Perhaps we should first define computational fluid dynamics. CFD is 
the numerical solution of a discrete approximation to the partial differential 
equations describing the fluid flow past a body. The governing equations may 
be potential, Euler, Reynolds averaged Navier-Stokes or boundary layer, de- 
pending upon the important flow phenomena. Usually finite difference or 
finite volume methods are used, but other discrete approximations (finite 
element, spectral) may be adopted. For subsonic or supersonic inviscid flow 
the potential equation can also be solved by panel methods. For an aircraft 
configuration these computations take minutes on the biggest computers 
available. However, the entire solution of the flow field is obtained, 
including the aerodynamic forces. 

Although the development of accurate, reliable, and economical CFD 
methods is far from complete, the ones which are available today can already 
treat many practical problems. Two-dimensional transonic airfoils are 
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routinely analyzed on the computer using methods pioneered by Bauer, 

2 

Garabedian, Korn and Jameson # Resulting airfoils are then tested in the 

wind tunnel to verify the design and to establish off-design performance. 

For 3-D configurations, panel methods are well developed and widely used for 

subsonic and supersonic analysis. The 3-D transonic methods are rapidly 

developing into reliable and accurate codes 3 . By way of illustration, F-l 

4 

shows calculations by Boppe and Aidala for the shuttle launch configuration, 

5 

F-2 presents calculations by Yu for the effect of a body on the transonic 
wing pressure distribution for a transport configuration, and F-3 shows the 
shock pattern cn a transonic compressor rotor calculated by Thompkins 6 . It 
can only be expected that in upcoming years the development of better algo- 
rithms and computers will produce improved methods for analyzing complex 
aerodynamic problems. 

The problem of interest to us here is the utilization of emerging CFD 
capability for aerodynamic design. A number of approaches and techniques 
are being developed or used. The easiest and most commonly found one is 
a numerical " cut- and- try" procedure using a direct analysis code. As new 
vehicle geometries are conceived, the aerodynamic performance is computed. 

In other words, CFD codes are used instead of wind tunnel tests to arrive 
at an acceptable design. This procedure has been incorporated in the design 
of several new aircraft with substantial savings of wind tunnel testing 
hours. F-4 illustrates this for the Grumman Gulf stream III. 

A second approach to the design problem is the "inverse method." As 
new aerodynamic understanding is developed, a pressure distribution is con- 
ceived and the geometry which corresponds to it is computed. Recent develop- 
ments of these techniques are given in R-7-11. One difficulty with inverse 
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techniques is insuring that realistic geometries are obtained for the 
specified pressure distribution. There are theoretical limitations for 
two-dimensional flow on what constraints may be specified 12 . Generally, 
the requirement of exactly matching the specified pressure distribution is 
replaced by a relaxed requirement which yields a practical geometry and 
some least squares best match to the specified pressure distribution. 

Another uncertainty in inverse methods is knowing which pressure distribution 
is really optimum. For example, should it contain a shock or be shockless? 
Garabedian and McFadden 11 found that specifying a shockless pressure dis- 
tribution on a swept wing did not lead to a shockless flow field solution. 
"Hanging shocks" were found to exist above the body which weakened before 
reaching the body. An optimization procedure was then used to modify the 
surface pressure distribution to achieve nearly shockless flow. 

There are a number of design techniques which are used to varying 
degrees but apply to a limited class of problems due to the basic assumptions 
inherent in the formulation. For example, hodograph methods 2,12 are extremely 
useful for two-dimensional, shock— free shapes but they are not extendable to 
more general problems. Chin and Rizzetta*s streamfunction method can solve 
two— dimens iona. inverse problems with shocks, but extension to three-dimensional 

problems is not clear. The clever artificial gas method of Sobieczky is 

applicable to two- and three-dimensional shock-free shapes 15 . Other methods 
exist which we have not mentioned. 

Another approach to design is to seek optimum shapes which satisfy a 
stated design objective. This approach was highly developed using calculus 
of variations for problems which could be treated by classical theoretical 
aerodynamics. The book by Miele represents the extent to which this approach 

has been developed. An "optimum body" is a shape which has the maximum or 



minimum value of a desired objective function, subject to specified constraint 

For example, elliptic wing loading yields the minimum induced drag for a given 

span, and a Sears-Haack axi symmetric body has the minimum wave drag for a 

fixed volume. The classical theory of optimum shapes often does not produce 

geometries which can be used directly in actual flight vehicles. However, 

they do serve as a standard of the best which could be achieved under ideal 

circumstances — a sort of Carnot efficiency of aerodynamics. 

A somewhat different approach to seeking optimum aerodynamic shapes 

which is suited to the CFD era is the use of numerical optimization methods. 

Generally speaking, numerical optimization procedures are search algorithms 

which seek to minimize (or maximize) an objective function subject to the 

specified constraints by systematically varying the free constants (design 

variables) which parameterize the system. Numerical optimization methods 

have been applied to many aerospace fields including structures, trajectories, 
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guidance and control, propulsion, and preliminary design . A recent article 
18 

by Vanderplaats gives a historical development of optimization methods in 

structures. The first application of numerical optimization to aerodynamic 

19 

problems is that of Hague, Rozendall and Woodward using linear panel 

methods. However, the most recent applications stem from the paper of 

20 

Hicks, Murman and Vanderplaats using transonic potential flow methods 

for airfoils. The farthest extensions of this work are the three-dimensional 

studies given in R-21-23. Also about the time the work of Hicks et al 

24 

appeared. Parsons, Goodson and Goldschmied applied the method to laminar 
flow control for underwater bodies. 

Our concern in this article is to examine the application of 
numerical optimization to aerodynamic design. The attractiveness of 
this approach lies in its general applicability* to practical "real world" 



configurations. The aerodynamic analysis can be as complicated as a 
reliable CFD code can handle. The constraints may include not only geometric 
restrictions, but flow properties and off design performance. On the other 
hand limitations arise from computational efficiency and the ability to 
describe the body by a few parameters. It is our contention in this article 
that Aerodynamic Design by Numerical Optimization (ADNO) is potentially the 
most powerful design method using CFD. We will describe the method in more 
detail and cite illustrative examples. However, it is clear that the method 
is not yet practical for routine applications and we will conclude with 
suggestions for key areas requiring research. 

In order to illustrate the numerical optimization technique, let us 
consider a hypothetical isolated wing design problem. First we define an 
objective function OBJ which is a single number representing the quantity 
we want to minimize. An example of an OBJ is the sum of the wave, friction, 
and induced drag. The definition of the objective function is up to the 
designer, but it must be expressible as a single number. If the designer 
wants to maximize some objective function, then OBJ is selected as the 
negative of that function. A frequently used objective function has been 
to minimize the deviation from a specified surface pressure distribution 
on the airfoil or wing. In this way an inverse problem is recast as an 
optimization problem. 

Second we must define design variables which are free constants to 
be selected in order to minimize the objective function. The usual choices 
are coefficients of some functions describing the body shape. For example, 
if we start with a baseline wing, perturbations on different shapes could 
be added to the thickness and camber distributions at selected stations. 

The amplitudes of these perturbation shape functions are the design variables 
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The wing twist and taper can also be design variables. In principle as 
many design variables (M) as desired may be used, but in practice the com- 
putational time increases approximately as M 2 . Thus the selection of design 
variables becomes the challenge of finding as few free constants which 
describe as broad a class of body shapes as possible. 

The M— dimensional space of design variables is called the desig n space. 
At each point in design space OBJ has a value. One can visualize this space 
with a topographical analogy, the contours of constant OBJ being elevation 
lines. The space is filled with hills and valleys, ridges and passes, and 
maybe even cliffs and cornices. The goal is to find the lowest point in 
design space. For a complicated problem, one can expect many local or 
relative minima. Cliffs correspond to discontinuities and cornices to non- 
uniqueness, both of which are bound to give mathematical difficulties. 

There is, of course, no guarantee that the flow field solution is unique. 
Recently, non— uniqueness has been demonstrated for transonic potential flow 
past an airfoil 25 . Let us consider the two-dimensional design space shown 
schematically in F-5. The design variables are a representative thickness 
ratio and taper ratio, and the objective function is some arbitrary performance 
parameter . 

All designs have constraints which must be enforced. For our illus- 
tration, consider geometric constraints, aerodynamic constraints, and off- 
design constraints. These may be represented as barriers in design space 
as illustrated in F-6. Geometric constraints dictated by structures might 
correspond to minimum and maximum values of taper ratio and a minimum thick- 
ness. These are horizontal and vertical lines in design space or linear 
constraints. An aerodynamic constraint might be that the maximum adverse 
pressure gradient parameter on the wing be below a specified ceiling to 
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avoid separation and buffeting# Off design constraints might reguire no 
leading edge separation at take-off and that the wing root bending moment in 
a maneuver is limited. These latter constraints are shown as curved barriers 
in design space as they most likely wxll be nonlinear functions of the 
design variables. For each selection of the design variables, CFD codes or 
other techniques must be available to evaluate the constraints and the 
objective function. 

The constrained design space for our illustration is the region 
inside the constraint boundaries of F-6 and is called the feasible region. 

The optimization problem is to find the point of minimum objective function 
MIN within the feasible region. Numerical optimization procedures are 
automated search processes which start from some initial point and seek 
the minimum. Since the topology of the design space is not known ahead 
of time, this is a "searching in the dark" task. There are a number of 
algorithms available^' ^ including random walk, non-gradient, and gradient 
methods to aid in this task. A complete review is outside the scope of 
this paper. 

The most widely used optimization method for aerodynamic design has 

28 

been the CONMIN algorithm developed by Vanderplaats . It is a gradient 
type constrained minimization algorithm based upon the method of feasible 
directions. Starting at an initial point in design space, gradients of 
OBJ are calculated. If no constraints are being violated, a search is made 
in the direction of steepest descent until a new minimum is found or a con- 
straint encountered. If the starting point is near an active constraint, 
the direction normal to the constraint surface must be found so that the 
search can be made in the direction which minimizes OBJ and stays in the 
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feasible region. At the conclusion of the search step, a new point in design 
space is established, new gradients are calculated, and another search per- 
formed, this time using a conjugate gradient method in place of steepest 
descent. The process continues until no reduction in OBJ is found. A more 
detailed description is given in R-29. 

There are several general comments which can be made regarding the 
numerical optimization approach. First, the calculations can be lengthy 
v/ith many searches and concomitant objective and constraint evaluations. 
Usually only a few constraints are "active" and that reduces the burden 
somewhat. Nevertheless, computational efficiencies of the algorithms are 
of critical importance. Second, the results are no more reliable than the 
aerodynamic models. Thus, design by optimization is no good unless the CFD 
codes have been thoroughly validated. Third, there is no guarantee a global 
minimum will be found. Although theoretically this may be a perplexing 
problem, it is of little practical importance. Each result is usually an 
improvement, and different starting points may be tried to seek new minimums. 
Finally, the method is very general. With the exception of the present man- 
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in-loop procedures, we are not aware of any other method which can incorporate 
off-design constraints. 

A number of publications have reported numerical optimization studies 
for low speed, • multi-element, transonic and circulation controlled airfoils, 
transonic wings, a propeller, and a low drag underwater body. We will cite 
only representative results to illustrate various features of the method. 

In the majority of the studies, the optimization and CFD analysis codes were 
coupled as illustrated in F-7. The optimization code is the "driver." It 
passes values of the design variables to the CFD analysis code. In turn the 
CFD code returns values of the objective function and aerodynamic constraints. 
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If gradients are required, they are calculated from one-sided finite 

difference formulas formed by individually perturbing each design variable 

about the center point. The optimizer then evaluates gradients and search 

directions, giving new values of design variables to the CFD code, etc. 

For M design variables, (M + l) CFD calculations are needed to compute 

gradients. Order M searches are required, each taking about three CFD cal- 

2 

culations. Thus, order 3M CFD evaluations are needed. Usually an iterative 
method is used in the CFD calculation to solve the flow equations. Since most 
of the aerodynamic evaluations represent small changes from a previous solution, 
convergence is much faster than for a regular analysis calculation 20 . Thus, 
the total computational time is much less than 3 M 2 regular solutions. 

Hicks, Human, and Vanderplaats considered the problem of minimizing 
shockwave drag at transonic speeds. The CFD code was based on inviscid 
transonic small disturbance theory for symmetric non-lifting bodies. Although 
this is not a practical design problem, some interesting features were illus- 
trated. F-8 shows a result where a virtually shock-free airfoil was obtained 
starting from an arbitrary airfoil with a sizeable shockwave. Off design 
drag results indicate a rapid increase past the design point. When an 
additional constraint dictating off— design behavior was enforced (drag rise 
between M= .8 and M = .81 less than 0.0015), the optimum airfoil has a weak 
shock (F-9) . F-10 illustrates an example for a supersonic flow with a bow 

shockwave. In this case, the minimum drag body has a detached shockwave due 

_ to the blunt nose. It is interesting that this body has a slightly lower 
drag than a parabolic arc of revolution with the same enclosed volume, which 
linear theory predicts should be the minimum drag body 16 . 

As a second example, we describe a study involving minimizing viscous 
drag • The application was to axisymmetric low drag under-water bodies 
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enclosing a specified volume. Drag minimization was obtained by delaying 

transition and avoiding separation. In this case the CFD model was a coupled 

potential flow panel method, boundary layer analysis, and standard transition 

prediction procedure. A polynomial description of the body with nine design 

variables was used and various geometric constraints were enforced. The 

constrained minimization procedure was of the random search variety. Results 

of the design indicated as "Body X-35" are shown in F-ll where they are 

compared to well known low drag body called the "Dolphin." The optimization 

procedure produced a body with a drag coefficient approximately 30% lower 

7 

than the Dolphin for a Reynolds number of 10 (F-lla) . Post-design calcu- 
lations indicate this favorable effect is realized over a wide Reynolds 
number range (F-llb) . 

The basic procedure of the Hicks, Murman, Vanderplaats study was 

29 

extended to more practical inviscid transonic flow calculations by using 

a CFD program based on the transonic full potential equation for lifting 

30 

airfoils. Researchers at Lockheed-Georgia Co. applied the method to a 
redesign study of the C141 airfoil. Two redesigns were done. The first 
involved modifying only the first 12% of the upper surface. The second 
involved modifying the entire upper surface. In both studies, the objective 
was to reduce the drag at the cruise Mach number of 0.74. For the second 
modification the design was constrained to have the same maximum thickness 
and nose bluntness as the baseline C141 airfoil, approximately the same lift 
coefficient, and also an off-design constraint that the wave drag at M=.72 
be less than 0.0020. This was to avoid "drag-creep" which can result from 
doing a single point design. 

For the first study where the leading edge region was modified, 
viscous effects are negligible and the design methodology was borne out by 
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wind tunnel verification (F-12) . The researchers 
test data showed that a 7% improvement in (ML/d) 


state: "The wind tunnel 
may have resulted from 


the modification of the upper surface of the forward 12% of the airfoil. 
In addition to producing an efficient airfoil modification, numerical 
optimisation required about half the computational time and resulted in a 


25% reduction in engineering hours when compared to a conventional trial- 
and-correction approach." For the second study involving modification of 
the entire upper surface, viscous effects are important. An attempt was 
made to correct for them but the result was unsuccessful. The wind tunnel 
tests showed a higher drag than the base line airfoil (F-12) . The negative 
result of the second study was due to the inability of the CFD coda to model 


the viscous effects. In other studies (e.g. R-22) where a simple displace- 
ment thickness was added, the airfoil design process was in better agreemen 
with experimental data. 


Several three-dimensional numerical optimization wing designs have 

2i_2 3 

been reported . The study by Haney, Johnson, and Hicks 22 sought to 
redesign the A-7 wing to produce substantial improvements. The CFD code 
used was an inviscid transonic full potential equation method (FL022) with 
a simple two-dimensional viscous displacement surface added. OBJ for this 
case was the deviation from a specified surface pressure distribution; that 
is, an inverse problem was solved. A total of 120 design variables were 
used. At each of five span stations, the twist, trailing edge camber, and 
twenty-two surface shape functions were used as design variables. It was 
prohibitive to vary all of these design variables simultaneously and instead, 
each section was optimized sequentially starting at the root and moving to 
the tip. The design goal of the first study (wing No. 1) was to increase 
the thickness by 71% with the same drag divergence Mach number. The design 
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goal of the second study (Wing No. 2) was to increase the thickness by 28%, 
reduce the induced drag by 25% , and maintain the same drag divergence Mach 
number. These design goals together with two-dimensional airfoil tests led 
to the specification of the desired wing pressure distribution which formed 
the objective function. The planfom shapes were preset by preliminary 
design studies. F-13 shows the planform geometries and results of the wind 
tunnel tests. Both wings show substantial improvement in drag characteristics 
compared to the original A- 7 wing without any degradation in drag divergence 
Mach number. The results of the numerical optimization procedure were 
successful. Although one of the three-dimensional inverse methods 7 ” 11 may 
have been more computationally efficient, the application of numerical 
optimization was relatively straightforward for this case. 

As a final example, we cite the results of a U.S. Air Force sponsored 
23 

study with Lcckheed-Georgia Co. and Grumman Aerospace Corp. The purpose 
was to develop and validate a new transonic wing design procedure using 
numerical optimization techniques.” Both fighter and transport configurations 
were considered in a three-phase study involving design, test, and evaluation. 
Although many limitations were encountered, the concluding remarks from the 
transport study include: "We have developed a new transonic wing design 

method using the numerical optimization scheme. We have also shown that 
new computational methods offer a means for the aerodynamic design of wings 
with transonic performance superior to that which could be obtained using 
previous design techniques. The method is relatively easy to use, and it 
is compatible with established industry design procedures. By using the 



new method, a 40% to 50% reduction in the cost associated with wing cruise 
aerodynamics design is obtainable. " The latter statement was made despite an 
estimated 5-10 hour computation time on a CDC-7600 computer. The fighter 
configuration study resulted in the statement: "The results described 

here show both the benefits and hazards of numerical optimization for aero- 
dynamic design. The numerical optimization will work best when both the 
flowfield analysis and numerical model of tne design problem arc accurate. 
Although greater computer resources are generally needed for more complex 
and accurate. analysis, the cost of numerical optimization would still be 
less than that of additional wind tunnel testing. As both computer and 
analysis code capability increase, numerical optimization will take a 

greater role in aerodynamic design." 

The above examples illustrate many advantageous features of aerodynamic 
design by numerical optimization (ADNO) and call attention to some of the 
limitations. Clearly the basic CFD model must produce reliable results 
before optimization can be used. However, with the enormous effort being 
expended to develop such codes, it can be expected that the future is bright 
in this category. The efficiency of the procedure depends in part on using 
aerodynamic "smartness" to specify the objective function, constraint function, 
and design variables. It is desirable to pick functions which are aero- 
dynamically "well conditioned." Perhaps some sort of adaptive parameterization 
can be developed to select those design variables which are making the biggest 
contribution and "turn off" the unimportant ones. Another possibility might 
be to further exploit the power of inverse methods by combining them with 
optimization techniques as suggested by Garabedian 31 . 

The most important limitation at present appears to be the large com- 
puter resources which are required. However, we note the experience of the 
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Lockheed-Georgia 30 team that if a careful accounting is done of the total 
engineer and computer hours to do a design, ADNO may not be as expensive 
as it appears. With the yearly reductions in computer costs and increase 
in computer capability, this situation should become less severe. Be 
that as it may, it is desirable to develop algorithms which are more 
efficient. The cited A DMO calculations have been performed by taking 
existing CFD and optimization codes and combining them in "black-box" 
fashion (F-7) . By far, the most time is spent in the CFD code doing aero- 
dynamic calculations. Two suggestions can be made for more efficient 
approaches. 

The first is to eliminate the requirement to execute the complete 
CFD calculation for every aerodynamic evaluation and to replace it by a 
suitable approximate solution (F-14) . For gradient calculations and as 
the minimum is approached, each new set of design variables should produce 
only a small change in the aerodynamic solution. Thus the complete CFD 
code need only be used for calculations producing large changes and an 
approximate calculation may be done for small changes. The selection 

criteria may be preset or done adaptively* Some demonstrations along this 

32 .33 

line are reported by Vanderplaats , Stahara, Elliot, and Spreiter , 

34 35 

Bristow, Hawk, and Thomas , and Peeters . 

The second suggestion is to combine the optimization and CFD 
algorithms into a one-step procedure in which the only converged solution 
is the final answer. One technique for solving large sets of discrete 
equations is to minimize the residuals. By using Lagrange multipliers, 
objective functions and constraints might be incorporated directly into 
the residual minimization. Labrujere*^ reports results using this approach. 



It Is hoped thst this article will stimulate interest and creative 
thinking shout the enormous potential of optimisation procedures for eero. 
dynamic design and ho» it can be developed into a practical approach. 
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F-8. Wave drag minimization for transonic airfoil with constraints 
on minimum thickness and enclosed volume. (a) Initial and 
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F-12. Measured airfoil performance for baseline 0141 airfoil compared with 
two modifications designed by numerical optimization (R-30) 
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